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Abstract 

Upward-going stopping muons initiated by atmospheric and 77^ interactions in the rock 
below the Soudan 2 detector have been isolated, together with a companion sample of neutrino- 
induced single muons, created within the detector, which travel downwards and exit. The 
downward-going sample is consistent with the atmospheric-neutrino flux prediction, but the 
upward-going sample exhibits a sizeable depletion. Both are consistent with previously reported 
Soudan-2 neutrino-oscillation results. Inclusion of the two samples in an all-event likelihood anal- 
ysis, using recent 3D-atmospheric-neutrino-flux calculations, reduces both the allowed oscillation 
parameter region and the probability of the no-oscillation hypothesis. 



PACS numbers: 14.60.Lm, 14.60.Pq, 96.40.Tv 



1 Introduction 

The recently published neutrino oscillation analysis of the Soudan-2 atmospheric-neutrino 
data [1] used only Fully Contained (FC) events and Partially Contained (PC) events in 
which the neutrino interaction vertex was contained within the Soudan 2 tracking calorime- 
ter. Not utilized by that study were two additional, topologically very similar, categories 
for which further analysis was needed to separate the samples and to eliminate their non- 
neutrino background. Using nomenclature introduced by the MACRO experiment [2], the 
two categories are labelled as UpStop events and InDown events. 

UpStop events are upward-going stopping muons which arise from charged-current and 

interactions occurring in the rock surrounding the Soudan 2 cavern; only the final state 
muon is detected as a stopping, non-interacting track in the detector. The muon may be 
accompanied by a small number of hits from a decay electron near its stopping point. 

InDown events are and interactions in the detector yielding downward going, exiting 
muons with three or fewer hits arising from hadronic track(s) at the production vertex. 
Approximately 65% of these InDown events are quasi-elastic interactions with low energy 
protons. Interactions having more than three hadronic hits at their primary vertex were 
included in the PC sample analyzed previously [1]. 

Separation of these two neutrino event samples is made possible by the fine-grain imaging 
of the Soudan-2 honeycomb-lattice tracking calorimeter. Their angular distributions exhibit 
features which are indicative of atmospheric oscillations. The samples can be 

incorporated in a straightforward way into the likelihood analysis described in Rcf. [1]. Their 
inclusion has enabled an improved determination of the — oscillation parameters from 
this experiment and a more stringent rejection of the no oscillation hypothesis. 

Analyses of upward through-going muons [2, 3, 4, 5], and stopping muons [2, 4, 6, 7] 
have been previously reported. Whereas through-going muon samples originate from a broad 
high-energy neutrino spectrum having a mean Ey of approximately 100 GeV, UpStop events 
originate predominantly from interactions with 1 20 GeV. Consequently they provide 

different constraints for oscillation scenarios. Among the underground tracking calorimeter 
experiments, MACRO has provided the most detailed treatment to date of UpStop and 
InDown events. In that experiment it was not possible to separate the two categories, so they 
were analyzed as a combined sample. Clearly, it is advantageous to separate the samples, 
since comparison of their zenith angle distributions can provide additional discrimination 
between low (rs 10~^ eV^) versus high (> 10""^ eV^) values of Am^. 

2 Detector and data exposure 

Soudan 2 was a 963 metric ton (770 tons fiducial) iron tracking calorimeter with a honeycomb 
geometry which operated as a time-projection chamber. The detector was located at a depth 
of 2070 meters-water-equivalent on the 27th level of the Soudan Underground Mine State 
Park in northern Minnesota. The calorimeter started data taking in April 1989 and ceased 
operation in June 2001 by which time a total exposure (fiducial exposure) of 7.36 kton-years 
(5.90 kton-years) had been obtained. 
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The calorimeter's tracking elements were 1 m long, 1.5 cm diameter hytrel plastic drift 
tubes filled with an argon-C02 gas mixture. The tubes were encased in a honeycomb matrix 
of 1.6 mm thick corrugated steel plates. Electrons deposited in the gas by the passage of 
charged particles drifted to the tube ends under the influence of an electric field. At the 
tube ends the electrons were amplified by vertical anode wires which read out a column 
of tubes. A horizontal cathode strip read out the induced charge and the third coordinate 
was provided by the drift time. The ionization deposited was measured by the anode pulse 
height. The calorimeter produced three-dimensional track hits with a spatial resolution of 
approximately 1 cm^ and separated by an average of about 3 mm of steel. The corrugated 
plates, interleaved with drift tubes, were stacked to form 1x1x2.5 m^, 4.3-ton modules from 
which the calorimeter was assembled in building- block fashion [8]. 

Surrounding the tracking calorimeter on all sides, but mounted on the cavern surfaces 
and well separated from the outer surfaces of the calorimeter, was a 1700 m^ Veto Shield 
array of two or three layers of proportional tubes [9]. The shield recorded the presence of 
cosmic ray muons coincident in time with events in the main calorimeter and thus identified 
background events, either produced directly by the muons or initiated by secondary particles 
coming from muon interactions in the rock walls of the cavern. Additionally, for neutrino- 
induced muons which enter or exit the tracking calorimeter, the shield array recorded the 
muon in time-coincidence with the event in the central detector. 

3 Separation of UpStop and InDown Events 

The event imaging afforded by the Soudan 2 tracking calorimeter made it possible to dis- 
tinguish the topologically similar UpStop and InDown events. Events of both types were 
isolated during routine processing. 

Events were classified as UpStop candidates if they satisfied the following criteria: 

(1) The track was muon-like, devoid of kinks or scatter vertices. 

(2) The track length was greater than 100 cm. 

(3) The muon endpoint occured in a live detector region. An event was removed if the 
candidate endpoint occured in the inactive region between modules. 

(4) Track ionization and straggling were consistent with the hypothesis of an upward-going 
muon which ranges to stopping. That is, near the edge of the detector the track was 
straight and lightly ionizing while, near the interior end, the track exhibited multiple 
scattering and/or heavy ionization. 

(5) Associated hits at the track endpoint, if any, had to be consistent with an electron 
shower from muon decay. 

An anode-versus-cathode view of an UpStop data event is shown in Fig. 1 where multiple 
scattering can be discerned as the muon ranges to stopping. Endpoint decay hits, the three 
hits modestly displaced from the muon endpoint in Fig. 1, are observed in some events (with 
higher probability for than for /i" since the former do not undergo nuclear absorption 
within iron nuclei). 
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Figure 1: An UpStop data event recorded in the anode-cathode matched view (front view of the calorimeter). 
Typically, multiple scattering becomes pronounced as the muon approaches its range endpoint. The endpoint 
decay shower of three hits favors identification as a /x+rather than a /U~ . 



InDown events were required to satisfy criteria 1-3 above and in addition: 

(4) The muon track was straight and hghtly ionizing at its interior end. Near the edge of 
the detector the track might, but need not, have exhibited ranging behavior in the form 
of muhiple scattering and/or heavy ionization. 

(5) Associated hits near the interior end of the track, if present, must have been consistent 
with hits from a proton or tt"^ track, lying in a straight hnc and heavily ionizing. 

These topological features are exhibited by the InDown data event shown in Fig. 2. At 
the event vertex, the muon is accompanied by a track of two hits for which the ionization is 
relatively heavy. This pattern is typical of a recoil proton. The event of Fig. 2 is a candidate 
quasi-elastic z/^ n ^ fJ'~P- 

There were a few events which could not be resolved as UpStop or InDown, the direction 
of the track being undetermined. Fortunately, ambiguous cases were rare for tracks which 
have visible lengths exceeding one meter in the detector. For the purpose of analysis, such 
events were retained as an Ambiguous category. 

4 Event processing and simulation 

Both UpStop and InDown data events were selected in the Partially Contained event sample, 
and were processed as described in Ref. [1]. The Monte Carlo sample of the contained- vertex 
InDown events was also part of the routine data processing, in which Monte-Carlo events 



4 



Soudon2 Data 

Run 1 12498 Event 1454 
07-Dec-2000 00:51:15.94 



195 1 



B 

o 



870 



Z (cm) 




1070 



Figure 2: An InDown data event recorded in the cathode versus time view (calorimeter side view). The 

muon emerges from an event vertex which is well contained; its trajectory, a straight line initially, undergoes 
small angle deflection as the ranging muon approaches the detector floor. A proton recoil of two hits is 
visible at the vertex. 



were inserted into and processed together with the data stream, their identity only being 
revealed in the final analysis stage. However, additional simulations, not included in the 
main data processing, were needed for the neutrino interactions in the rock surrounding the 
cavern, which give rise to UpStop events [10]. 

4.1 Simulation of upwEird-stopping muon events 

The GEANT Monte Carlo program together with modified Soudan-2 software provided the 
UpStop simulation (UpStop-MC). A total of 68.7 million neutrino interactions in greenstone 
rock were simulated. The event vertices were distributed randomly through rock volumes 
which were centered on the Soudan 2 cavern. Since high-energy charged-current (CC) events 
can project muons to the cavern from more remote rock than low-energy events, the dimen- 
sions of the primary rock volumes were chosen to increase with increasing E^,. Final-state 
particles were tracked through the rock by GEANT and fourvectors of particles that reached 
the veto-shield array were saved. These were then passed through the Soudan-2 Monte Garlo 
to produce realistic detector hits superimposed on detector noise represented by random- 
trigger records. 

UpStop-MC events were then processed through the standard Soudan-2 triggering, recon- 
struction and selection software for PG events. For UpStop-MG events yielding ionization 
within the tracking calorimeter (71,000 events), the survival rates decreased monotonically 
with increasing primary Ei,, refiecting the diminishing probability for muons from energetic 
events to stop in the detector. Survival rates ranged from 12.8% for events with E,, < 10 
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GeV, to 2.5% for events with > 40 GeV. 

4.2 UpStop cuts and scanning 

The PC selection filter required candidate tracks to penetrate to the fiducial region while 
not being through-going. Additional requirements, detailed in Ref. [1], were imposed in 
order to reject the high-flux background of downward-going cosmic-ray muons. A total of 
7662 UpStop-MC events passed the filter and simulated trigger requirements. However, only 
34% of these events yielded a potentially interesting topology in the detector. Consequently, 
additional cuts were applied to the true kinematic variables to reject those events which 
were certain not to pass the subsequent analysis cuts. These cuts (existence of a final-state 
muon with cosine zenith angle, cos 9^ < +0.05, and energy upon arrival at the detector 
within the range 350 MeV < E^ < 3500 MeV) were designed to ensure that the event had 
an upward-going muon that stopped within the calorimeter fiducial volume [10]. A 54% 
sample of the surviving events was then scanned by physicists, using scanning rules identical 
to those used for PC data event scanning. The additional criteria given in Sec. 3 were also 
applied to both data and MC events. All events which satisfied the scanning criteria were 
then reconstructed manually using the experiment's standard interactive graphics software. 
Only the reconstructed sample was used in the subsequent analysis. It contained a factor of 
25 more events than the data sample. 

4.3 MC event rate normalization 

The atmospheric neutrino flux used to generate the UpStop events was the onc-dimcnsional 
calculation of the Bartol group [11], modulated by the solar cycle as described in Ref. [1]. 
Other fluxes were simulated by applying correctional weights to the generated events. For 
consistency with Ref. [1], the numbers and plots in Sects. 5 and 6 were weighted to cor- 
respond to the updated ID Bartol-96 flux [12]. The oscillation analysis of Sect. 7 of this 
paper used the latest three-dimensional fluxes from the Bartol group [13] and Battistoni et 
al. [14]. The neutrino cross sections were those encoded in NEUGEN3 [15]. The target 
nuclear composition was that of Soudan rock, described in Ref. [16]. The effect of Pauli 
blocking in elastic and quasi-elastic reactions was accounted for, however nuclear effects on 
resonance production and on deep inelastic scattering flnal states were neglected. 

The event rate calculations have a sizable systematic error. For the comparison of this 
data with the MC presented in Sections 5 and 6, a normalization factor of 0.85, determined 
from the measured rate, assuming no oscillations, was applied. In the oscillation analysis 
described in Sec. 7 the overall normalization was a free "nuisance" parameter. 

5 Event rates and backgrounds 

5.1 Backgrounds in UpStop events 

Two sources of background events were considered: 
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(1) Cosmic-ray muons which scatter in the rock and eventually enter the detector in an 
upward direction. 



(2) Charged hadronic tracks, especially pions, produced at large angles in interactions of 
cosmic-ray muons in the rock surrounding the detector. 

Unlike experiments which are situated under mountains, the flat overburden at Soudan 
ensures that the flux of cosmic-ray muons becomes less than the flux of neutrino-produced 
muons signiflcantly above horizontal angles [16]. Thus the background from the flrst source 
is negligible. 

Hadronic interactions of neutrons produced in cosmic ray interactions were shown in Ref. 
[1] to be a background to contained neutrino events. A similar, related, flux of charged 
hadrons also emerges from the cavern walls and enters the detector. 

There are two distinguishing features of entering hadronic tracks: 

(1) Veto-Shield signal: In addition to the Veto Shield hit corresponding to the passage of 
the hadron track, there arc likely to be extra hits due to other particles produced in 
the muon interaction. In general, the track in the detector will not be aligned with 
these extra hits. It is thus useful to distinguish between the total number of in-time 
Veto Shield hits (nj^f ) and the number of in-time hits geometrically associated with the 
incoming/leaving track (n]^f). 

(2) Range: Hadronic tracks have a limited range in the detector before stopping or inter- 
acting. 
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Figure 3: Distributions of Veto-Shield hit counts versus track range, (a) for UpStop data candidates and (b) 
for UpStop Monte Carlo events, prior to muon length cuts. 

Fig. 3 shows n^f versus the range of the stopping track for each event of the UpStop data 
(Fig. 3a) and of the MC (Fig. 3b). A clear excess of data events with large n^f and small 
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Figure 4: Distribution of track range for candidate UpStop data events (crosses) compared to the neutrino 

Monte Carlo sample (shaded histogram), for events with visible track length exceeding one meter. The 
distributions (normalized to the events with > 500 g/cm^) are shown (a) before and (b) after the requirement 



vs 

''all 



VS 
Hrk 



has been applied. Fig. 4b shows good agreement between the data and the neutrino MC. 



range is observed, corresponding to the expectation for incoming hadronic tracks. However, 
at track lengths beyond the hadronic range, the data is consistent with the UpStop MC. The 
background signature is emphasized in Fig. 4a which shows the projection onto the range 
axis of the UpStop data (points with errors) and the MC (shaded histogram). The MC is 
normalized to the data with range > 500 g/cm^ (> 3.8 pion interaction lengths). Fig. 4b 



shows the same distributions but with the additional constraint that n 



vs 

all 



VS 
'hrk 1 



i.e. all 



Veto-Shield hits must be geometrically associated with the incoming track. The MC is then 
in good agreement with the data. However, to ensure that the residual hadronic background 
is negligible, a cut requiring the track range to be greater than 260 g/cm^, corresponding to 
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two pion interaction lengths, is also applied. Since the calorimeter is, to good approximation, 
a uniform medium of 1.6 g/cvo? density, the range requirement corresponds to a minimum 
track length requirement of ~160 cm. The effective muon momentum threshold for UpStop 
and InDown events is > 530 MeV/c. 

Finally, the cosine of the zenith angle, cos^^, of the reconstructed UpStop track is required 
to be smaller than +0.05. 



5.2 UpStop backgrounds using hadronic scatter events 

The PC data analysis also recorded events with an incoming track making a hadronic scatter. 
A sample of 25 data events was obtained which can be used to gauge the background from 
incoming non-scattering hadrons. A representative event is the upward-going, stopping, 
charged pion track shown in Fig. 5. There are two coincident hits in the shield which are 
in close proximity to the track's entrance point into the cavern, hence nj^f — 2. There is an 
additional coincident hit in the shield floor, hence n^f = 3. Inside the detector, the track 
scatters and stops. The scatter, which gives rise to a recoil proton, plus the track ionization, 
identifies the track as a tt^ rather than a /x^ or a proton. 
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Figure 5: A scattering upgoing pion track (data event) accompanied by three time-coincident Veto Shield 
hits. 

Distributions of cos^^ and of track range versus n^f for the incoming hadronic sample are 
shown in Fig. 6. Fig. 6a shows that most of the events are downgoing. Among the seven 
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Figure 6: Distributions of (a) cos6'2 and (b) the number of coincident veto shield hits versus the track range, 
for events which exhibit hadronic scatters. In Fig. 6(a), the shaded (open) histogram includes events where 
all the in-time Veto-Shield hits are (are not) associated with the detected track. 



upward-going tracks, five have n^if > n^^f. (open histogram); the remaining two events with 
n^if = nj^f (shaded histogram) are close to horizontal, and have short range. No hadronic 
events pass the shield and track length selections for upgoing tracks (ignoring their visible 
scatters). On the basis of this observation, plus the agreement in track range distributions 
between the UpStop data and the neutrino MC of Fig. 4, backgrounds from non-scattering 
charged pions and protons are estimated to contribute less than two events to the UpStop 
data and are hereafter neglected. 



5.3 Backgrounds in InDown events 

A potential background for InDown events arises from downward, through-going cosmic ray 
muons whose entrance into the detector was not detected due to a rare episode of poor or non- 
existent ionization drifting within a calorimeter module. Great care was taken to record all 
such incidents during data taking and additional checks were made by studying individual 
module efficiencies as a function of time. A special scan was carried out which rejected 
events for which there was a possibility of such an occurrence. Additional discrimination 
against this background was provided by the active shield array, since through-going muons 
yielded pairs of time-coincident hits having widely separated spatial locations. Using the 
same minimum track range as was used for the UpStop sample (> 260 g/cm^), no InDown 
events had n^f > nj^f . Other backgrounds for PC events were shown to be negligible in Ref. 
[1], thus the InDown sample was assumed to be background- free in the analysis presented 
below. 
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5.4 Event rates 



A final sample of 1081 fully reconstructed UpStop-MC events was retained for subsequent 
analysis. Within this simulation sample, 80% of events originate with neutrino interactions 
having E^, < 10 GeV. Charged-current quasi-elastic scattering accounts for about one third 
of the interactions. The neutrino fraction / {v^ + of the sample is 64%. 

The numbers of candidate data neutrino events which survive are listed in Table 1, where 
the MC numbers have been scaled to an exposure of 5.90 fiducial kton-years and include the 
factor of 0.85 to normalize to the event rate of Ref. [1]. Comparison of the data with 
the neutrino MC predictions, the sum of columns 2 and 3, shows that the observed InDown 
rate is consistent with the prediction, whereas the UpStop data rate appears suppressed by 
a factor of approximately two. These trends are in agreement with the expectation from the 
oscillation analysis of Ref. [1]. Note also the small size of the ambiguous sample and the 
small misidentification rate between the UpStop and InDown samples. 

Tabic 1: Numbers of data and Monte Carlo events which pass all cuts. The no-oscillation MC event rate is 
normalized to the measured c-flavor event rate of Ref. [1] assuming no oscillations. 



Assigned as 


No-osc. MC Truth 


Data 


InDown 


Ui)Stoj) 


InDown 


12.4±1.4 


0.3±0.1 


16 


UpStop 


1.8±0.5 


53.3±1.8 


26 


Ambig 


0.8±0.3 


3.4±0.4 


2 



6 Energy and angular distributions of UpStop/InDown neutrinos 

The neutrino energy, E,,, for Monte Carlo UpStop and InDown events is shown in Fig. 7. 
The UpStop events have an average E^ of 6.2 GeV. In contrast, the InDown events have 
lower El, values with an average of 2.4 GeV. 
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Figure 7: Comparison of the primary spectra for UpStop and InDown events. 



11 



The muon track provides a good estimator for the incident neutrino direction. For UpStop 
events the average angle between the incoming neutrino and the muon is 11°. For the lower 
energy InDown sample, the average angle is 13°. 

In the likelihood analysis of Ref. [1], the variable logiQ{L/E) was used, where L is the 
neutrino pathlength in kilometers and E is the neutrino final state visible energy in GeV. 
As demonstrated in Ref. [1], multiple scattering of the exiting muon in PC events allows 
sufficient neutrino energy determination to permit the reconstruction of the \ogio{L/E) dis- 
tribution. Fig. 8 shows the \ogiQ{L/E) distribution for InDown data (crosses) compared to 
the no-oscillation MC. The solid-line histogram represents the sum of the partially-contained 
InDown events from the contained- vertex PC Monte Carlo (PC-MC, dark shading), plus 
misidentified UpStop-MC events (hght shading). Good agreement is observed between the 
data and the unoscillated Monte Carlo. Hence no evidence is found for oscillations in 
atmospheric neutrinos which are incident from directions above the horizon. 
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Figure 8: Distribution of \ogin{L/E) for the InDown events. The data are shown as crosses; the soHd-line 
histogram shows the neutrino MC expectation for no oscillations. The MC prediction arises mostly from 
contained vertex PC events (PC-MC, dark shading), but receives a small contribution from UpStop-MC 
events which were misidentified as InDowns (light shading). 



For the UpStop events, however, the whole of the hadronic-shower energy plus a fraction of 
the muon energy is missing. Thus, whereas L can be calculated accurately from cos6z, E is es- 
sentially undetermined. The average difference between \ogio{L/ E)^, and logio{L/ E)i^ is 0.69 
(FWHM=0.66), spanning four bins used in the oscillation analysis. Therefore logio{L/E) 
is not a useful variable for analysis of oscillations in these events. In this section we show 
data as a function of cos6'^. In the analysis described in Section 7 it is convenient to fit 
as a function of logiol/. Fig. 9 shows the distribution of cos^^ for upward-stopping muon 
data events. The distribution decreases steadily towards the nadir. Fig. 9 also includes 
the simulated UpStop-MC (light shading) and the misidentified contained-vertex PC-MC 
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Figure 9: Distribution of cos6z for the UpStop events. The erosses represent the data, and the sohd-hne 
histogram shows the expected MC distribution for no oscillations. Light and dark-shaded areas show the 
contributions of UpStop and InDown events for no oscillations. 

(dark shading) distributions, for the no-oscillation case. Significant disagreement between 
the neutrino UpStop events and the no-oscillation expectation is apparent towards the nadir, 
which is consistent with the loss of upward-going /^.-flavor events due to oscillations. 

The distributions of Figs. 8 and 9 imply constraints on neutrino oscillation scenarios. 
These samples have been included in a likelihood analysis together with all the other neutrino 
events from the experiment. The method and the results of this global fit to Soudan-2 
neutrino data are described in Section 7 below. 

7 Oscillation analysis 
7.1 Outline of the method 

The oscillation analysis is a bin-free hkelihood analysis based on the prescription of Feldman 
and Cousins [17]. A detailed description of the method can be found in Ref. [1] and is not 
repeated here. The main points of the analysis are: 

(1) As reported in Ref. [7] and confirmed by Ref. [1], the distributions of neutrino-induced 
e-flavor data events are consistent with the null oscillation MC predictions, up to an 
overall normalization. Only the /x-fiavor data exhibit oscillation effects. Thus this 
analysis assumes two flavor z/^ Ur oscillations. 

(2) The FC and PC samples described in Ref. [1] are used unchanged in the present analysis. 
The InDown events are added to the PC muon-flavor sample, and the Ambiguous events 
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are used only in the overall normalization. The UpStop events are treated as a new 
category and analyzed as a function of logiol/ rather than as a function of \ogio{L/ E). 

(3) A likelihood function for the data is constructed as a function of Am^, the mass-squared 
difference, and of sin^ 29, where 9 is the mixing angle, using probability density functions 
(pdf's) determined from the MC sample. Details of the formahsm are given in Ref. [1]. 

(4) The summed negative log likehhood is evaluated at each point on a 15 x 80 grid of 

sin^ 2^^ X logio(Am^) with sin^ 2^ varied between 0.0 and 1.0 and Am^ varied between 
10^^ and 10*^ cV^. The lowest negative log likelihood on the grid is found and AC, the 
difference between the lowest value and the value in each (sin^ 29, Am^ ) grid square, is 
plotted. 

(5) A background contribution of non-neutrino events arising from neutrons and gammas 
produced by muon interactions in the rock around the detector is added to the hkelihood 
function. The background contribution only affects the FC events; the PC events and 
the new InDown and UpStop events are treated as background-free. 

(6) The overall normalization of the MC and the amounts of background in the different 

FC event samples, estimated using shield-tagged data events and the depth distribution 
of the event vertices, are nuisance parameters whose values are optimized at each grid 
square. 

(7) The allowed confidence level regions are calculated by the method of Feldman and 
Cousins [17]. That is, MC experiments are generated and analyzed at each grid square to 
calculate the expected hkelihood rise for a given confidence level based on the statistical 
and systematic errors at that grid square. In addition to the systematic errors described 
in Ref. [1], a 10% systematic error on the relative normalization of the UpStop events 
to the remainder of the data was allowed. The latter error represents uncertainties 
in density and nuclear composition of rock below the detector, and uncertainties with 
variation of neutrino cross sections in rock versus iron. 

(8) The analysis of Ref. [1] used the one-dimensional fiux calculation of the Bartol group 
[12]. This analysis uses their new three-dimensional calculation [13] and compares it 
with their ID calculation and with the 3D calculation of Battistoni et al. [14]. 

The 44 new data events documented here are added to the 488 data events of the previous 
analysis. However, the new events are "high resolution" /i-fiavor events, those most sensitive 
to oscillations, consequently they enhance the sensitivity afforded by the 167 events of that 
type in the previous analysis. 

7.2 Oscillation results 

The values of AC are plotted in Fig. 10 as a function of sin^ 26* and logio(Am^). The 
resulting surface is similar to that reported in Ref. [1]. The main difference is that the 
likelihood rise for the grid square with the lowest values of Arri^ and sin^ 29 (called the 
no oscillation point) is 16.0 compared with 11.3 for the previous analysis. The new data 
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Figure 10: The data likelihood difference, A£, plotted as a function of sin 29 and logio(ATO2). 

has significantly increased the discrimination against the no-oscillation hypothesis, mostly 
because of the large suppression of the UpStop data compared to the MC prediction. The 
probability of the validity of the no-oscillation hypothesis is discussed in Sect. 7.3. 

The A£ surface in Fig. 10 exhibits two nearly equal minima, one at the grid square 
centered at Am^ = 0.0017 eV^, sin^ 2^ = 0.97, and one at Am^ = 0.0052 eV^, sin^ 2^ 
= 0.97. The first minimum is the "best fit point" of this analysis, while the second minimum 
was the best fit point in the previous analysis [1]. The difference of A£ between the two is 
only 0.18. There is a small rise in likelihood between the two minima which peaks at about 
1.8 in the region of the Super-K best-fit point. However, the 90%-confidence-level limit of 
this analysis, determined in Sect. 7.3, contains all of the Super-K allowed region. The value 
of the flux normalization at the best fit point is 91% of the Bartol-3D prediction [13]. 

Although the analysis is carried out using the log-likelihood function, it is useful to eval- 
uate the goodness of the fit by projecting out the distributions for the various data sets and 
calculating a for the data compared to the MC prediction. Fig. 11a shows the data for 
the total muon PC sample (including the InDown events) and Fig. lib shows the UpStop 
events, compared to the MC predictions. 

The PC events are in good agreement with the no-oscillation histogram, but disagree, 
particularly at low L/E (downward going u events), with the prediction at the highest 
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Figure 11: (a) The logj^Q L/E distribution for the PC event sample that includes the InDown events, (b) The 
UpStop events plotted as a function of logioi- The points with error bars are the data. The solid histogram 
is the MC prediction at the best fit point, and the dashed histogram shows the no-oscillation expectation. 
The dotted histogram depicts "saturated oscillations" with sin^ 29 = 1.0, Am^ =1.0 eV^. 



considered values of Am^ and sin^ 26* (~1.0 eV^ and ~1.0), where the oscillations have the 
greatest frequency and largest size (called "saturated oscillations"). On the other hand, the 
UpStop events are in disagreement with the no-oscillation prediction, particularly at large L, 
but agree with the saturated oscillation prediction. This behavior is indicative of oscillations 
with values of Am^ in the low 10~^ eV^ region where downward-going neutrinos have not yet 
oscillated and upward-going neutrinos have saturated oscillations. The best-fit prediction 
agrees well with both distributions. 

Table 2 gives the fo^' the comparison of data and MC for the PC and UpStop distribu- 
tions as well as the full data set including the other distributions described in Ref. [1]. The 

at the best fit point, using all of the data, is 35.6 for 30 data points. 



7.3 Confidence levels 

The 90%-confidence-level surface ACgo (not shown), generated by the MC experiments under 
Feldman- Cousins prescription is very similar to that of the previous analysis. Combining the 
ACgo surface with the data likelihood surface of Fig. 10 gives the solid-line 90%-confidence- 
level contour shown in Fig. 12a. For comparison, the 90%-confidence-level contour deter- 
mined previously in Ref. [1] is shown by the dotted line. It can be seen that the allowed 
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Table 2: The for comparisons of the data to various MC predictions for the PC events, the UpStop events, 
and the full data set. 





X^/Number of data points 


PC 


UpStop 


All data 


Best fit 

No oscillations 
Saturated oscillations 


5.0/5 
5.9/5 
19.8/5 


0.6/4 
7.9/4 
1.9/4 


35.6/30 
66.6/30 
63.2/30 




0.4 0.6 0.8 1 0.4 0.6 0.8 1 

sin2(2e) sin2(20) 

Figure 12: (a) The 90% confidence allowed region from the Feldman-Cousins analysis of this work (solid 
line), compared to that of Ref. [1] [1] (dotted Hue), (b) Contours at 68%, 90%, and 95% confidence level 
(dotted line, thick solid line, and dashed line respectively), compared to the 90% sensitivity contour (thin 
solid line). 
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region is more restrictive and that some of the contour structure indicated by the previous 
analysis has been smoothed. This is due to the fact that the hkehhood surface is rather fiat 
at the base of the valley and small changes in the data can move the contours substantially 
in this region. 

Fig. 12b shows the 68%, 90% and 95% contours, a further indication of the shape of the 
hkelihood surface. Also displayed in Fig. 12b is the 90% sensitivity contour obtained from the 
Monte Carlo experiments, which denotes the expected 90% contour for experiments with this 
statistical precision and systematic errors. As was the case in Ref. [1], the 90%-CL contour 
from this analysis is more restrictive than the estimated sensitivity contour, due, in part, to 
a small mismatch of the overall event-rate normalization in the electron and muon samples. 

The probabihty of no oscillations can be calculated, under the Feldman-Cousins prescrip- 
tion and including all of the systematic effects, by generating experiments at the minimum 
Am^ and sin^ 26 grid square and counting those MC experiments that give a larger likelihood 
difference than 16.0. In 500,000 simulated experiments, 16 had a larger likelihood difference 
giving a probability of 3.2 x 10~^ for the no-oscillation hypothesis. 

For UpStop muons and upward-going contained vertex events, the initiating neutrino 
may experience Mikheev-Smirnov-Wolfenstein (MSW) resonance and other matter effects 
as a result of the traversal of thousands of kilometers of terrestrial matter. The magnitude 
of matter-induced deviations from vacuum oscillations was studied using simulated event 
samples; the samples were weighted in accordance with three-neutrino mixing and the nor- 
mal mass hierarchy, using the approximation of a uniform (path-weighted mean) terrestrial 
density [19]. For the range of plausible Am^ values, it was found that matter effects, even 
with maximally allowed mixing angles, can only introduce a few percent additional depletion 
of muon-flavor neutrinos beyond that which results from i/^ — > v^. vacuum oscillations [20] . 
Since the scale of such deviations is well below the statistical sensitivity afforded by the data, 
matter effects were neglected in this analysis. 

7.4 Flux model comparison and event rate normalization 

The analysis of this data has been carried out for three different atmospheric neutrino fiux 
calculations: (i) the one-dimensional fiux calculation of the Bartol group [12], (ii) the Bartol 
three-dimensional calculation [13], and (in) the three-dimensional calculation of Battistoni 
et al. [14]. As well as the flux prediction, the analysis requires an estimate of the height 
in the atmosphere at which the neutrino is produced. This is particularly important for 
neutrinos coming from overhead where the path in the atmosphere is a large fraction of 
the total path length. A parameterization of the pion and muon decay heights was made 
using the formalism of Ref. [18] for the one-dimensional Bartol calculation. A similar 
parameterization for decay heights was prepared for the three-dimensional Bartol case. The 
latter parameterization was also used for predictions based upon the Battistoni et al. flux. 

Fig. 13 shows the 90%-CL region for the three cases. There is a small change from the 
one- to three-dimensional flux models, however the three-dimensional models of Bartol and 
Battistoni et al. give almost identical results. The only significant difference between the 
three cases is in the absolute normalization of the flux. At the best-flt point of this analysis, 
the normalization factor (number of events observed/calculated) for the Bartol ID flux is 
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Figure 13: The 90%-confidence-level allowed region for three atmospheric-neutrino flux calculations. The 
dotted curve is based upon the ID-model of the Bartol group from 1996 [12]. The more recent 3D calculations 
of Battistoni et al. [14] and of the Bartol group [13] lead to the dashed and solid curve, respectively. 

0.86, while for the Bartol 3D flux the factor is 0.91 and for the Battistoni et al. flux it is 



The authors of the flux calculations typically quote large errors of ±20% on the absolute 
normalization, due to the uncertainties on the incoming cosmic-ray fluxes and on nucleus-air 
cross sections. There are also significant errors on the neutrino cross sections. It is thus of 
interest to determine the experimental error on the ratio of the measured to the predicted 
event rate. The experimental event rate is proportional to the incident neutrino flux, the 
neutrino cross sections in the detector, and the detector acceptance. This experiment can 
thus determine the normalization of the atmospheric neutrino flux at the Soudan-2 site times 
the neutrino cross sections encoded in the NEUGEN3 program [15], for an iron calorimeter 
with a given energy threshold. Translation of this normalization factor to other experiments 
at the Soudan site and at other sites is possible in principle. However it requires knowledge 
of the relative neutrino fluxes at the different sites and the ratio of the neutrino cross sections 
if a different detector medium or a different neutrino generator is used. 

The event-rate normalization factor for Soudan 2 is subject to the following errors: 

(1) The total number of neutrino events observed in this experiment above an energy thresh- 
old of 300 MeV, obtained from Table 1 of Ref. [1] and Table 2 of this paper, is 481.2±26.2 
(±5.4%). The error includes the statistical error and the error on the background sub- 
traction. 

(2) The statistical error on the Monte Carlo sample is ±1.6%. 

(3) The variation of the fitted normalization factor over the 68% confidence region of the 
oscillation parameters is ±3%. The Feldman-Cousins analysis includes the systematic 



1.02. 
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errors associated with the background subtraction, cross section uncertainties and en- 
ergy scale uncertainties. 

(4) Depending on the value of 62^-, Ve ~^ oscillations with the solar parameters could 
change the flux of Uf, that have traversed the Earth [21, 22]. The change can be positive 
or negative depending on whether 623 is smaller or greater than 45°. Using the Super- 
K limits for sin^ 2^^23 and recent values for the solar oscillation parameters [7, 23], an 
uncertainty of ±3.3% in the calculated electron-event rate is estimated. 

(5) Uncertainty arises in the rate of multi-GcV muon events due to matter effects [20]. A 
±1.6% error contribution to the event rate calculation is inferred. 

(6) Any mismatch between the Monte Carlo representation of the detector and reality 
could introduce a relative error in the acceptance of the two and thus an error in the 
normalization ratio. Detailed studies of individual event channels revealed no signiflcant 

differences [24, 25]. The relative proportions of different event types and event rejection 
modes in data versus Monte Carlo samples were studied. A ±2% systematic error, 
estimated from the maximum differences found between the data and Monte Carlo, has 
been assigned to account for uncertainties arising from geometric acceptance and other 
residual effects. 

Based upon the ±8% quadrature sum of these errors, an overall normalization factor of 
0.91±0.07 is determined for this analysis. This normalization is specific to the Soudan 2 
site, the detector medium, the Bartol 3D flux, and to the neutrino cross sections encoded in 
the NEUGEN3 event generator. It is averaged over the years 1989-2001, one full solar cycle. 
The same percentage error, ±8%, is apphcable to the Soudan-2 normalization factors given 
above for the Bartol ID and Battistoni et al. 3D atmospheric fluxes. 

As a check, a normalization which is mostly independent of the oscillation pa- 

rameters can be obtained from the total electron sample of contained and partially contained 
events. Table I of Ref. [1] lists 208.7±15.9 background-subtracted electron-neutrino events 
to be compared with an expected rate from the 3D Bartol prediction of 238.1 events, yielding 
a normalization factor of 0.88±0.07, where the error is just statistical from the number of 
events and does not include the contribution from v^. oscillations or the other error 

sources detailed above. 

8 Conclusions 

Samples of upward stopping muons produced by neutrino interactions in the rock below 
the Soudan-2 detector and partially-contained events with downward-going muons produced 
in the detector have been separately isolated. These two new data sets provide additional 
support and constraints for the hypothesis of atmospheric-neutrino oscillations. The flux of 
upward stopping neutrino-induced muon events is observed to be suppressed by a factor of 
approximately two, while downward-going muon events are not suppressed. An oscillation 
analysis using the method described in Ref. [1] and adding this new data gives a more 
restrictive 90%-confldence-level allowed region of Am^ and sin^ 29 . The probabihty of the 
no-oscillation hypothesis is reduced by more than a factor of 10, to 3.2 x 10~^. 
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Figure 14: The Soudan 2 90% confidence allowed region in sin^ 29, l\m? (solid line) compared with the 
allowed regions of MACRO (dot-and-dashed line) [27], and of the Super-K zenith angle [7] (dotted line) and 
L/E [26] (dashed line) analyses. 



The data have been analyzed using three models of the atmospheric flux at the northern 
geomagnetic latitude of this experiment. The models include two recent three-dimensional 
flux calculations and an older one-dimensional calculation. The oscillation parameters are 
found to be essentially independent of the flux calculation. The normalization factor for 
the experiment, 0.91 ±0.0 7, is the measured event rate divided by the calculated event rate 
where the latter is the convolution of neutrino fluxes of the Bartol 3D flux calculation with 
neutrino cross sections encoded in NEUGEN3. The denominator for this ratio contains 
elements which are specific to the Soudan-2 detector analysis, and site. Consequently, the 
normalization factor cannot be compared in a straightforward way to other experiments 
at other geomagnetic latitudes with different detector media and using different neutrino 
interaction generators. 

Comparison of this experiment's revised 90% CL allowed region with the most recent 
Super-K [7, 26] and MACRO [27] allowed regions is shown in Fig. 14. This result is in good 
agreement with both experiments. 
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